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Abstract. Seismic electric signals have been found to obey the ubiquitous l/f a 
behavior [Phys. Rev. E 66, 011902(2002)]. The newly introduced concept of natural 
time enables the study of the dynamic evolution of a complex system and identifies 
when the system enters the critical stage. On the basis of this concept, a simple model 
is proposed here which exhibits the 1//° behavior with a close to unity. Furthermore, 
we present recent data of electric signals, which when analyzed in the natural time 
domain are found to exhibit critical dynamics and hence can be classified as seismic 
electric signals. 
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1. Introduction 

Among the different features that characterize complex physical systems, the most 
ubiquitous is the presence of l/f a noise in fluctuating physical variablespQ. This 
means that the Fourier power spectrum S(f) of fluctuations scales with frequency 
/ as S(f) ~ 1//°. The power-law behavior often persists over several orders of 
magnitude with cutoffs present at both high and low frequencies. Typical values 
of the exponent a approximately range between 0.8 and 4 (e.g., see Ref.[2] and 
references therein), but in a loose terminology all these systems are said to exhibit 
1/f "noise". Such a "noise" is found in a large variety of systems, e.g., condensed 
matter systems(e.g. [3]), freeway traffic [U El [6], granular flow [7], DNA sequence [H], 
heartbeat [9] , ionic current fluctuations in membrane channels 1 10], river discharge [TT], 
the number of stocks traded daily p2], chaotic quantum systems [T31 [141 fT5l fT6] . the light 
of quasarsJTTJ, human cognition[18] and coordination [19], burst errors in communication 
systems [20 J, electrical measurements [2T], the electric noise in carbon nanotubes[22j and 
in nanoparticle films [2 3], the occurrence of earthquakes |24j etc. In some of these systems, 
the exponent a was reported to be very close to 1, but good quality data supporting 
such a value exist in a few of them[3], e.g., the voltage fluctuations when current flows 
through a resistor [25]. 

The l/f a behavior has been well understood on the basis of dynamic scaling 
observed at equilibrium critical points where the power-law correlations in time stem 
from the infinite-range correlations in space (see Ref.[2] and references therein). Most 
of the observations mentioned above, however, refer to nonequilibrium phenomena for 
which -despite some challenging theoretical attempts [261 1271 [281 129] - possible generic 
mechanisms leading to scale invariant fluctuations have not yet been identified. In other 
words, despite its ubiquity, there is no yet universal explanation about the phenomenon 
of the 1/ f a behavior. Opinions have been expressed (e.g., see Ref.[T3]) that it does not 
arise as a consequence of particular physical interactions, but it is a generic manifestation 
of complex systems. 

It has been recently shown [301 131~1 [321 [331 EU ESI [36] that novel dynamic features 
hidden behind the time series of complex systems can emerge if we analyze them in terms 
of a newly introduced time domain, termed natural time x- It seems that this analysis 
enables the study of the dynamic evolution of a complex system and identifies when 
the system enters a critical stage. Natural time domain is optimal [37] for enhancing 
the signal's localization in the time frequency space, which conforms to the desire to 
reduce uncertainty and extract signal information as much as possible. The scope of 
the present paper is twofold. First, a simple model is proposed (Section 2) which, in 
the frame of natural time analysis, leads to 1/ f a behavior with an exponent close to 
unity. Second, we present the most recent experimental data on Seismic Electric Signals 
(SES) which are transient low frequency (< 1Hz) signals observed before earthquakes 
[351 1391 l40l [4lT [42] . since they are emitted when the stress in the focal region reaches 
a critical value before the failure [43J. The analysis of the original time series of the 
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SES activities have been shown to obey a 1/ /-behavior [31]. Here in Section 3, we show 
that the SES activities observed in Greece during the last months exhibit the features 
suggested, on the basis of natural time, to describe critical dynamics. Such features 
have been found[H] for laboratory data [4"5l 146] of the avalanches in a three dimensional 
pile of rice (which is similar to the prototype example of sandpiles used in the proposal 
of the Self-Organized Criticality, SOCj26][27]) getting progressively closer to the critical 
state. 

In a time series comprising N events, the natural time Xk = k/N serves as an 
index [30] I3TI [32] for the occurrence of the fc-th event. The evolution of the pair (xk, Qk) 
is studied [Sni EDE21I331E11E51E51I1ZISSII131ED!, where Q k denotes a quantity 
proportional to the energy released in the fc-th event. For example, for dichotomous 
signals, which is frequently the case of SES activities, Qk stands for the duration of the 
k-th pulse. The normalized power spectrum II (a?) = |$(u;)| 2 was introduced [30] 131"] . 
where 

( k\ 



k=l 

^N 



and Pk = Qk/ ^2 n= i Qn, u> = 27T0; <p stands for the natural frequency. The continuous 
function $(o;) should not be confused with the usual discrete Fourier transform, which 
considers only its values at <fi = 0, 1, 2, . . .. In natural time analysis, the properties of 
U.(u) or 11(0) are studied ([30] [31]) for natural frequencies less than 0.5, since in this 
range of 4>, II (a>) or 11(0) reduces ([301 SB E21 HE]) to a characteristic function for the 
probability distribution pk in the context of probability theory. When the system enters 
the critical stage, the following relation holds[30l EH 09]: 
„, , 18 6cosu; 12 sin uj 

n (oj ) = — — . 2 



For u -> 0, Eq.© leads to [30 



U(u) w 1 - 0.07u 2 



which reflects [49] that the variance of x is given by 
«i = (X 2 ) - (X) 2 = 0.07, 

where (f(x)) — Ylk=iPkf{Xk)- T^ 16 entropy S in the natural time-domain is defined 
as [IB [31] 

s = (x^x) - (x) ln (x), 

which depends on the sequential order of events [35] [36] and for infinitely ranged temporal 
correlations its value is smaller [34] 148] than the value S u (= In 2/2 — 1/4 w 0.0966) of 
a "uniform" distribution (as defined in Refs. [30l [33] Ell EH] EE], e.g. when all pk 
are equal or Qk are positive indepedent and identically distributed random variable of 
finite variance), i.e., S < S u . The same holds for the value of the entropy obtained [47] 
upon considering the time reversal T, i.e., Tpk = pN-k+i, which is labelled by S_. In 
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Figure 1. (color) (a):Example of the evolution of e„ (see the text) versus the number 
of renewals n, i.e., in natural time, (b): The Fourier power spectrum of (a); the 
(green) solid line corresponds to 1/f and was drawn as a guide to the eye. (c): The 
DFA of (a) that exhibits an exponent a, o fa very close to unity, as expected from 
(b). (d):Properties of the distribution of e„ that explain the small deviation of (Xbfa 
from unity. The average value (e n ) (plus) and the variance ((e — (e n )) ) (crosses) as 
a function of n. The straight solid line depicts ln(n) and was drawn for the sake of 
reader's convenience. 



summary, the SES activities, when analyzed in natural time exhibit infinitely ranged 
temporal correlations and obey the conditions: 

Ki = 0.07 (3) 

and 

S,S_<S U . (4) 

The validity of these conditions for the most recent electric field data will be investigated 
in Section 3. 

2. The model proposed 



Here, we present a simple competitive evolution model which results, when analyzed in 
natural time, to 1/ f a "noise" with a very close to unity. Let us consider the cardinality 
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Figure 2. (a)The six(=3!) equally probable outcomes after the selection of 3 random 
numbers by the same PDF. Actually, the sample space is (in one to one correspondence 
to) the permutations of 3 objects, (b) The 24(=4!) equally probable outcomes after the 
selection of 4 random numbers by the same PDF. Again, the sample space is (in one to 
one correspondence to) the permutations of 4 objects. For the readers convenience, in 
each outcome, the corresponding e„-value (n — 3 or 4)is written. An inspection of (b), 
shows that p(e 4 = 1) = l/4,p(e 4 = 2) = 11/24, p(e 4 = 3) = 1/4 and p(e 4 = 4) = 1/24. 



e n of the family of sets S n of successive extrema obtained from a given probability 
distribution function (PDF); S equals to the empty set. Each S n is obtained by following 
the procedure described below for n times. Select a random number rj n from a given 
PDF and compare it with all the numbers of S n _i- In order to construct the set S n , we 
disregard from the set S n _i all its members that are smaller than r] n and furthermore 
include r] n . Thus, S n is a finite set of real numbers whose members are always larger or 
equal to r) n . Moreover min[S n ] > min[S n _i] and max[S n ] > max[S n _i]. The cardinality 
e n = |S n | of these sets, which may be considered as equivalent to the dimensionality of 
the thresholds distribution in the coherent noise model (e.g. see Ref . [51] and references 
therein), if considered as time-series with respect to the natural number n (see Figfl^a), 
which was drawn by means of the exponential PDF) exhibits 1/ f a noise with a very close 
to unity, see Figfj](b). This very simple model whose evolution is depicted in Figflja), 
leads to a detrended fluctuation analysis[9] (DFA) exponent a^FA ~ 1-02, see Fig{T](c), 
being compatible with the 1/f power spectrum depicted in FigHJ^b). The mathematical 
model described above corresponds to an asymptotically non-stationary process, since 
(e n ) oc Inn with a variance ((e — (e n )) 2 ) oc Inn (see FigJIJd)), and this conforms to the 
fact that the DFA exponent is slightly larger than unity. 

We now discuss an analytical procedure which clarifies some properties of the model. 
In order to find analytically the distribution of the probabilities p(e n ), one has simply to 
consider the possible outcomes when drawing n random numbers r\ n . Since the selection 
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is made by a means of a PDF, all these numbers are different from each other, thus 
-when sorted they- are equivalent to n points (sites) lying on the real axis. The value of 
t n varies as {r) n } permutate along these n sites independently from the PDF used in the 
calculation. Thus, a detailed study of the permutation group of n objects can lead to an 
exact solution of the model. It is well known, however, that the number of the elements 
of this group is n! and this explains why we preferred to use the numerical calculation 
shown in FigJTJ Some exact results obtained by this method are the following: (ei) = 1; 
(e 2 ) = 1 + 1/2, since p(e 2 = 1) = p(e 2 = 2) = 1/2; (e 3 ) = 1 + 1/2 + 1/3, since 
p{e 3 = 1) = l/3,p(e 3 = 2) = 1/2 and p{e 3 = 3) = 1/6; (e 4 ) = 1 + 1/2 + 1/3 + 1/4 (see 
Figj2]). Figure [2] analyzes the results for n = 3 (Figj2|a)) and n = 4 (FigJ2(b)). One can 
see that the probability p(e n = m) equals to the sum of the n possible outcomes as i] n 
moves from the left to right in the n columns of Figj2[ In each column, the probability 
to have at the end e n = m is just equal to the probability to keep m — 1 numbers from 
the numbers already drawn that are larger than r] n . This results in 



Using Eq.flSJ), one can prove that (e n ) = (e n -i) + which reflects that (e n ) = 



3. The recent electric field data 

We now proceed to the presentation of the most recent experimental data of SES 
activities recorded in Greece by means of the procedure described in detail in 
Refs. [52| 1531 14*8] . Figure [3] depicts four electrical disturbances that have been recently 
recorded at three measuring stations termed Lamia (LAM), located at 150km north- 
west of Athens, Patras (PAT) « 160km west of Athens, and Keratea (KER) « 30km 
east south-east of Athens. The signals are presented here in normalized units, i.e., 
subtracting the mean value and dividing the result by the standard deviation. The two 
electric signals in Figj3](a) were recorded at LAM on November 13, 2006 and January 
8, 2007, respectively, the latter being of larger actual amplitude than the former. In 
the upper panel of Figj3]Jb) an electric signal recorded at PAT on February 8, 2007 is 
depicted, whereas in the lower panel we also insert -for the sake of comparison- the signal 
that was recorded at the same station almost one year ago ,i.e., on February 13, 2006. 
This recent signal, i.e., the one on February 8, 2007, has an amplitude « 70% larger [54"] 
than that recorded one year ago which has been analyzed in Ref . |55j . Its subsequent 
seismic activity in discussed in the Appendix A (while more recent SES activities are 
presented in Appendix B). Finally, in Figj3]^c) we depict a signal recorded at KER on 
February 11, 2007. 

All these four recent signals were analyzed in natural time and found to be consistent 
with the conditions (j3J) and (j4j), thus they can be classified [50] [55] as SES activities. 
For example, in the lowest panel of FigJH we depict the evolution of the parameters 
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Figure 3. Electric signals recorded at LAM (a), PAT (b) and KER (c). (sampling 
rate f e x P =l sample/sec). The actual electric field amplitude E (in mV/km) in (c) is 
appreciably smaller than that in (a) and (b). 



«i, S and S- for one of these signals, i.e., the one recorded at KER on February 11, 
2007. The leftmost panel of this figure shows the original time series (along with its 
dichotomous presentation), while the rightmost panel indicates how the signal is read in 
natural time. Note, in FigJU(c), that at x — 1 (i- e -> upon the completion of the signal) 
K\ reaches the value K\ = 0.068 ± 0.003, while S and >S_ are S = 0.071 ± 0.003 and 
S- = 0.078 ± 0.003. Thus, we have assured that, within the experimental uncertainty, 
K\ ~ 0.07 and S, S_ < S u . Furthermore, upon shuffling the Qk randomly, we have 
found that, in all these four signals, the variance k\ and the entropies S, S- turn to be 
equal to the values expected from a "uniform" distribution, which assures [55] that their 
self-similarity solely stems from temporal correlations. 



4. Conclusions 



In summary, using the newly introduced concept of natural time: (a) a simple model 
is proposed that exhibits 1/ f a behavior with a close to unity and (b) electric signals, 
recorded during the last few months in Greece, are classified as SES activities since 
they exhibit infinitely ranged temporal correlations. This excludes any possibility 
of attributing these signals to nearby man-made sources, because the latter have 
weaker temporal correlations (i.e., their Hurst exponent lies usually in the range 0.5 
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Figure 4. (color) The electric signal recorded at KER on February 11, 2007 (red) 
along with its dichotomous representation which is marked by the dotted (blue) line, 
(b) How the signal in (a) is read in natural time, (c) The values of «i, S and 5_ as 
the signals evolves from \ — to \ — 1 versus the natural time x- Note that upon 
the completion of the signal, i.e., at \ = 1, the K\ value is close to 0.07, whereas 
both S and S- are smaller than S u . For the reader's convenience, the red and blue 
horizontal lines show the values n u = 1/12 and S u = 0.0966 of K\ and S 1 , respectively 
that correspond to a "uniform" (it) distribution. 

to 0.75 [331 EH)- 

Appendix A. What happened after the SES activity at PAT on February 8, 
2007, until April 23, 2007. 

According to the Athens observatory (the seismic data of which will be used here), 
a series of strong earthquakes (EQs) with magnitudes ranging from 5.0 to 6.0-units 
occurred as follows: First, a 6.0 EQ at Kefallonia area, i.e., 38.34°N 20.42°E, at 13:57 
UT on March 25, 2007. Second, a cluster of four magnitude class 5.0 EQs on April 10, 
2007 with an epicenter close to Trichonida lake, i.e., around 38.5°N 21.6°E. 

We show that the occurrence time of the impending strong EQ activity can be 
estimated by following the procedure described in Refs. [30j HHJ |49j [55l [56] , as it was 
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indicated in Ref . [5"4] . (We clarify that, during the last decade, preseismic information 
based on SES activities is issued only when the magnitude of the strongest EQ of the 
impending EQ activity is estimated to be -by means of the SES amplitude[38-40, 52,53]- 
comparable to 6.0 units or larger [15] .) 

We study how the seismicity evolved after the recording of the SES activity at 
PAT on February 8, 2007, by considering either the area A-.N^qE^'q or its smaller area 
B : -^37.'6^20.'o> which surround the EQ epicenters and the PAT station (see Fig JAl( a)). 
If we set the natural time for seismicity zero at the initiation of the SES activity on 
February 8, 2007, we form time series of seismic events in natural time for various time 
windows as the number N of consecutive (small) EQs increases. We then compute the 
normalized power spectrum in natural time 11(0) for each of the time windows. Excerpt 
of these results, which refers to the values deduced during the period from 20:53:59 UT 
on March 19 to 11:56:30 UT on 25 March, 2007, is depicted in red in FigfA2](a). This 
figure corresponds to the area B with magnitude threshold (hereafter referring to the 
local magnitude ML or the 'duration' magnitude MD) M t hres = 3.2. In the same figure, 
we plot in blue the power spectrum obeying the relation (2) which holds, as mentioned, 
when the system enters the critical stage (u = 27T0, where stands for the natural 
frequency). The date and the time of the occurrence of each small earthquake (with 
magnitude exceeding (or equal to) the aforementioned threshold) that occurred in area 
B, is also written in red in each panel. An inspection of this figure reveals that the red 
line approaches the blue line as N increases and a coincidence occurs at the last small 
event which had a magnitude 3.2 and occurred at 11:56:30 UT on March 25, 2007, i.e., 
just two hours before the strong 6.0 EQ. To ensure that this coincidence is a true one 
(see also below) we also calculate the evolution of the quantities K\,S and S- and the 
results are depicted in Fig. IA2( b) and lA2( c) for the same magnitude thresholds for each 
of the areas B and A, respectively. 

The conditions for a coincidence to be considered as true are the following (e.g., see 
Refs. [30, 48, 49, 55,56]): First, the 'average' distance (D) between the empirical and 
the theoretical 11(0) (i.e., the red and the blue line, respectively, in Fig JA2( a)) should 
be smaller than 10~ 2 . See Fig. IA2( b).(c) where we plot (D) versus the conventional 
time for the aforementioned two areas B and A, respectively. Second, in the examples 
observed to date, a few events before the coincidence leading to the strong EQ, the 
evolving 11(0) has been found to approach that of the relation (2), i.e., the blue one in 
Fig JA2( a) , from below (cf. this reflects that during this approach the Ki-value decreases 
as the number of events increases). In addition, both values S and S_ should be smaller 
than S u at the coincidence. Finally, since the process concerned is self-similar (critical 
dynamics), the time of the occurrence of the (true) coincidence should not change, in 
principle, upon changing the (surrounding) area used in the calculation. Note that in 
Fig. IA2( b). upon the occurrence of the aforementioned last small event at 11:56:30 UT 
of March 25, 2007, in area B the (D) value becomes smaller than 10 -2 . The same was 
found to hold for the area A (see Fig JA2( c)). 

An inspection of Fig. IA2( c) shows that upon the occurrence of the four magnitude 
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Figure Al. The maps in (a) and (b) show the areas A,B and C,D, respectively. 
The stars in A stand for the strong 6.0 EQ that occurred on March 25, 2007 in 
Kefallonia (leftmost star) and the four magnitude 5.0-class EQs on April 10, 2007 
close to Trichonida lake. 



class 5.0 EQs on April 10, 2007, the parameters «i, S, and 5_ jumped to values larger 
than those of the "uniform" distribution and then started to decrease again. Among 
these parameters, however, S still remains appreciably higher than S u , while S_ it 
already became smaller than S u and K\ gradually approaches 0.070. 

In order to elucidate the situation, we proceeded to the additional study of two areas 
C and D, which did not include the epicenters of the EQs at Kefallonia and Trichonida 
lake, i.e., C: N^^E^'g and its smaller region D: Nffi^E^J. The corresponding 
parameters are depicted in Fig JA3l which exhibit the following interesting feature: They 
gradually decrease, thus approaching the aforementioned values that characterize the 
proximity to the critical point (CP). The occurrence of the next few small events in 
these areas will reveal whether CP will be finally reached or not. This study is still in 
progress. 

The interesting feature seen in Fig. IA3I seems to conform with the following fact: 
Just on April 23, 2007, a single SES (in the sense discussed in Ref. 52) was recorded 
(Fig. IA4( a)) which was followed by an SES activity (Fig. IA4( b)). Comparing the actual 
amplitude of this single SES to that of the SES activity recorded on February 8, 2007 
(Fig. IA4( c)). we find that it is evidently stronger. 

Appendix B. Further experimental data added on May 8,2007. 

At 01:34 UT on May 7,2007 an EQ of magnitude 5.1 occurred with an epicenter at 
37.62°N21.08°E depicted with a star in Fig lBl( a). It actually lies within the areas D,C 
(Fig JAl( b)) which exhibited the interesting feature discussed in advance in Fig JA31 

Concerning the electrical data: Beyond the electrical recordings on April 23, 2007 
depicted in Fig JA4( a).(b). two additional electric signals have been recorded at PAT on 
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Figure A2. (color) (a) The normalized power spectrum(red) 11(0) of the seismicity 
as it evolves event by event (whose date and time (UT) of occurrence are written in 
each panel) after the initiation of the SES activity on February 8, 2007. The excerpt 
presented here refers to the period 19 to 25 March, 2007 and corresponds to the area 
B Mthres = 3.2. In each case only the spectrum in the area cf) £ [0, 0.5] is depicted 
(separated by the vertical dotted lines), whereas the H((j>) of Eq.(2) is depicted by blue 
color. The minor horizontal ticks for (f> are marked every 0.1. (b), (c) Evolution of 
the parameters (D), k\, S and S- after the initiation of the SES activity on February 
8, 2007 for the areas B (Mthres — 3.2) and A(Mthres — 3.2), respectively. In (b) the 
period just before the 6.0 EQ is shown, whereas in (c) the period is extended until 
03:26:33.4 UT on April 23, 2007. 



April 24, 2007 (FigE^Ka)) and April 18, 2007 (FigjB^b)); which have an amplitude 
markedly larger than that on February 13, 2006 (see also Ref . |54j ) . These have been 
classified as SES activities after analyzing them in natural time and applying the relevant 
criteria. For example, if we read in natural time the signal on April 24, 2007 (Fig JB2( d)) - 
the dichotomous representation of which is marked by the dotted (blue) line in Fig JB2( c)- 
we find the values «i = 0.067 ± 0.003, S = 0.072 ± 0.003, S_ = 0.069 ± 0.003 which do 
obey the conditions (3) and (4). 

We now investigate the seismicity after the aforementioned three SES activities on 
April 18, April 23 and April 24, 2007. The investigation is made in the areas D and E 
(Fig JBl( b)). the latter being almost equal to the former region C but slightly extended 
to the east, i.e., E: Np^EjoJ. Starting the computation of seismicity from the initiation 
of the SES activity on April 18, 2007, we obtain the results depicted in Figs lB3( a) and 
IB3( b) for the areas D and E, respectively. An inspection of the parameters (£>), Ki,S and 
S- reveals that they exhibited a true coincidence (as discussed in Appendix A) around 



Natural time and 1/f "noise 



12 




Figure A3, (color) (a) and (b) depict the evolution of the parameters (D), K\, S and 
SL after the initiation of the SES activity on February 8, 2007, but for the areas C 
and D, respectively, until 00:49:12 UT on April 23, 2007. 
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Figure A4. (a): Electrical recordings at PAT on April 23, 2007. (b): Excerpt of (a), 
in an expanded time scale, to show the SES activity indicated by the arrow in (a), (c): 
For the sake of comparison, the SES activity recorded by the same station on February 
8, 2007, is depicted. 



April 30, 2007, i.e., around one week before the 5.1 EQ on May 7, 2007 mentioned 
above. This study of the evolution of seismicity in the areas D, E still continues. 

In view of the fact that three SES activities (April 18, 23 and 24, 2007) have 
been observed with different polarities (see Fig JA4( a) and Figs JB2( a).(b)) while only 
one strong EQ has occurred to date, our investigation is also currently extended to the 
area F: Ng^Egf'^ (which is a sub- region of E to the east of PAT). This study is still 
in progress in order to discriminate whether this area may approach CP or not. 

Appendix C. Further data until May 14,2007. 

As mentioned in the previous Appendix, we also currently study the area F:Npg 5 Egf'fl 
shown in Fig lCll Starting the calculation from the initiation of the SES on April 18, 
2007 (FigjB2|b)) the resulting parameters (until 02:22:17 UT on May 14, 2007) (D), k x , 
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20°E 21 "E 22°E 23°E 24°E 20°E 21 "E 22°E 23°E 24°E 



Figure Bl. The maps in (a) and (b) show the areas C,D and D,E, respectively. The 
map in (a) is a reproduction of Fig lAlf b) (originally submitted on April 24, 2007) 
which also includes the 5.1 EQ (star) on May 6, 2007, close to Zakynthos Island. In 
(b) the area F: Njjf^Ef^fs (not shown) is also currently studied. 
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Figure B2. (color) (a),(b): Electrical recordings at PAT on April 24, 2007, and April 
18, 2007 respectively, (c): The electric signal depicted in (a) in normalized units along 
with its dichotomous representation which is marked by the dotted (blue) line, (d): 
How the signal in (c) is read in natural time. 
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Figure B3. (color) (a) and (b) depict the evolution of the parameters (D), k%, S and 
S- after the initiation of the SES activity on April 18, 2007, but for the areas D and 
E, respectively, until 04:14:12 UT on May 6, 2007. 
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Figure CI. A map showing the areas G and F. 
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Figure C2. (color) (a) and (b) depict the evolution of the parameters (D), ki, S and 
5_ after the initiation of the SES activity on April 18, 2007, but for the areas F and 
G, respectively, for M thres > 3.0 until 02:22:17 UT on May 14, 2007. 
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S, and >S_ from the computation of seismicity are shown in Fig lC2( a). for M t h res > 3.0. 

To investigate the spatial invariance, we also give in Fig JC2( b) the relevant results for 

a subregion G of F, i.e., G:Np;f E|f :f 5 . Interestingly, both Figs. IC2Ta). (b) exhibit a 

systematic tendency to gradually approach CP. The extent to which they will finally 

obey the aforementioned conditions for a true coincidence or not, can be judged upon 
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